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ABSTRACT 

The Mercury-Man ganese star 4> Her is a well known spectroscopic binary that has been the subject 
of a recent study by I Zavala et al . (2006), in which they resolved the companion using long-baseline 
interferometry. The total mass of the binary is now fairly well established, but the combination of the 
spectroscopy with the astrometry has not resulted in individual masses consistent with the spectral 
types of the components. The motion of the center of light of 4> Her was clearly detected by the 
Hipparcos satellite. Here we make use of the Hipparcos intermediate data ('abscissa residuals') and 
show that by combining them in an optimal fashion with the interferometry the individual masses 
can be obtained reliably using only astrometry. We re-examine and then incorporate existing radial- 
velocity measurements into the orbital solution, obtaining improved masses of 3.05 ± 0.24 M a and 
1.614 ± 0.066 M Q that are consistent with the theoretical mass-luminosity relation from recent stellar 
evolution models. These mass determinations provide important information for the understanding 
of the nature of this peculiar class of stars. 

Subject headings: binaries: general — binaries: spectroscopic — methods: data analysis — stars: 
chemically peculiar — stars: fundamental parameters — stars: individual (4> Her) 



1. INTRODUCTION 



In a recent paper I Zavala et alj (|2006f ) reported inter- 
ferometric observations of the Mercury-Manganese star 
4> Her (HD 145389, HR 6023, HIP 79101, a J2 ooo = 
16 h 08 m 46?18, (5j 2 ooo = +44°56'05'/7, spectral type 
B9:p(HgMn), V = 4.22), which is a binary star. This 
object belongs to a class of peculiar non-magnetic B-type 
stars that show abundance anomalies of several elements, 
some of which (such as Hg and Mn) c an be enhanced 
by orders of magnitude (jPrestonl Il974h . Depletions of 
other elements are seen as well. These anomalies are 
thought to be produced b y radiatively-dr iven diffusion 
and gravi tational settling dMichaudl [l970h . The obser- 
vations of Za vala et al.l (|2006[ ) spatially resolve the com- 
panion of <j> Her for the first time. The object had previ- 
ously been known as a single-lin ed spectroscop i c binary 
with a period of about 560 days (|Babcocklll971t lAikmanl 
Il976| ) and an eccentric orbit. Their interferometric mea- 
surements with the Navy Prototype Op tical Interferom- 
eter (NPOI) allowed IZavala et alj (|20Q6f ) to measure the 
brightness of the companion and use it to estimate its 
spectral type. This, in turn, facilitated their detection 
of spectral lines of the secondary for the first time in 
spectra taken at the Dominion Astrophysical Observa- 
tor y. The secondary was also d etected spectroscopically 
bv lDworetskv fc Willatti pOOtl) . 

The system is therefore tech nically now d o uble- lined, 
although as pointed out by I Zavala et alj (|2006f ) the 
handful of secondary radial velocities they were able 
to measure with great difficulty do not provide a firm 
constraint on the velocity semiamplitude of that star. 
Furthermore, they found that the combination of their 
high-precision interferometric orbit with the elements 
of the s ingle- lined spectroscopic orbit reported by 
lAikmanl (1976) led to absolute masses for the corn- 
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ponents that are inconsistent with the spectral types. 
Therefore, although the total mass of the binary is 
now fairly well known, the individual masses cannot 
yet be determined dynamically. Such basic properties 
of the stars are of considerable interest given the 
chemical peculiarities of this class of objects, which 
have been the subject of ex tensive studies of many dif- 
feren t kind s (see, e.g.. lAdelman. Gulliver fc Ravld 

20011: lAdelman. Adclma nfe Pintado! 120031: 

Adelman et all 120041: iDolk. Wa hkrcn & Hubrigl 120031 : 



Dworetskv fc Budai [ 120001 : lLeushinNl995L and references 
therein) . 

We were puzzled by this apparent inconsistency be- 
tween the seemingly precise radial velocities and the as- 
trometry, and we wondered whether other velocity mea- 
surements available in the literature might clarify the 
situation. In addition, (f> Her was obse rved in the course 
of the Hipparcos mission (lESAlfT99l . Those measure- 
ments clearly revealed the motion of the center of light 
due to the binary orbit, suggesting they could be used 
to good advantage in the determination of the individual 
masses. The motivation for this paper is therefore three- 
fold: i) To show how, in the absence of spectroscopy, the 
astrometric measurements from Hipparcos can indeed be 
combi ned with the NPOI observations of [Zavala et al.l 
(2006) to yield reliable masses for both stars for the first 
time, purely astrometrically. This serves as an interest- 
ing example of the value of the Hipparcos intermediate 
data for solving or improving binary orbits and inferring 
other stellar p roperties; ii) To readdress the issue of the 
lAikmanl (|1976f ) velocities in light of the new solution. We 
show that those data are not really inconsistent with the 
astrometry, but do lack critical phase coverage and must 
be supplemented by other information in order to be use- 
ful; Hi) To incorporate additional velocity measurements 
not previously used, in order to further strengthen the 
orbital solution. The final masses not only have much 
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improved precision, but are consistent with the mass- 
luminosity relation as given by current stellar evolution 
models. 

2. INTERFEROMETRIC OBSERVATIONS 

The origi nal NPOI measurements obtained by 

IZavala et al.l (2006) consist of interferometric squared 
visibilities (V 2 ) and closure phases collected on 25 sepa- 
rate nights from 1997 April to 2005 July. The visibilities 
of each night were combined to infer the angular separa- 
tion (p) and position angle (9) of the binary, and these 
were subsequently used to derive the elements of the as- 
trometric orbit. The measurements published are p and 
9, along with the corresponding error ellipses. The resid- 
uals from the orbit are only a few tenths of a milli-arc 
second (mas) in angular separation, and a few tenths 
of a degree in position angle, indicating the very high 
precision of those measurements. The orbit has an an- 
gular semimajor axis of 32.1 mas, which, when combined 
with th e Hipparcos paralla x (7THi P = 14.27 ± 0.52 mas) 
allowed [Zavala et alJ (|2006l ) to infer the total mass of the 
binary as M to t — 4.7 ± 0.6 M Q . Additionally, the inter- 
ferometric measurements yielded the average magnitude 
difference between the primary and secondary of cj> Her 
at two wavelengths: Am (5500 A) = 2.57 ±0.05 mag and 
Am (7000 A) = 2.39±0.05 mag. As discussed below this 
information is key to the mass determinations. 

3. HIPPARCOS OBSERVATIONS 

The star (f> Her was observed by the Hipparcos satel- 
lite during its 37-month astrometric mission under the 
designation HIP 79101, and was measured a total of 76 
times from 1989 December to 1993 March, or a slightly 
over two orbital cycles of the binary. Each measurement 
consisted of a one-dimensional position ('abscissa') along 
a great circle representing the scanning direction of the 
satellite, tied to an absolute frame of reference known 
as the International Celestial Reference System (ICRS). 
In addition to providing the trigonometric parallax as 
well as the position and proper motion of the center of 
light based on the synthesis of all these measurements 
(the so-called standard 5-parameter fit) , the astrom et- 
ric solution reported in the Catalogue ( ESA] [19971 ) re- 
vealed a perturbation that was modeled as orbital mo- 
tion, since the binary nature of the object was already 
known. Several of the orbital elements including the pe- 
riod (P), the eccentricity (e), and the longitude of peri- 
astron of the primary (u>a) were held fixed at t he values 
determined spectroscopically by I Aikmanl (|1976l ) to facil- 
itate the solution, and the other elements were solved 
for. These included the semimajor axis of the center of 
light (ciphot = 9.09 ± 0.65 mas), the inclination angle 
(i = 36° ± 14°), the position angle of the ascending node 
(fi = 188° ± 12°, J2000) and the time of per iastron pas- 
sage (T). More recently I Jancart et al.l (|2005l ) reported a 
reanalysis of the Hipparcos intermediate data, still adopt- 
ing P, e, loai a s well as T from the spectroscopic work by 
Aikman (1976), but with the added assumption that the 
secondary contributes no light to the system. This pro- 
vides a connection between the measured semimajor axis 
Qphot and the radial v elocity semiamp litude Ka, which 
was adopted also from I Aikmanl (|1976[ L Their results for 
the inclination (i = 10?3±0?7) and especially the position 
angle of the node (fi = 148?3 ± 3?0) are quite different 



from the Hipparcos solution, b ut we now kn o w fro m the 
detection of the secondary by I Zavala et al.l |2006) that 
th is star does contribu te some light, so the assumption 
of I Jancart et al.l (|2005l ) is not a good one in this partic- 
ular case. 

As it turns out, the information provided by the 
Hipparcos data is complementary to that given by the 
ground-based interferometry, and as we show below the 
combination of the two allows the individual masses of 
4> Her to be determined independently of any assump- 
tions, and independently also of any spectroscopic in- 
for mation. Since some of the orbital elements reported 
by IZavala et al . ( 2006) are much improved compared to 
lAikmanl (|1976f ). a reanalysis of the Hipparcos data with 
those constraints would therefore seem to be in order. 
Better still, the Hipparcos data can be combined directly 
with the NPOI observations in a simultaneous least- 
squares solution yielding all the elements at once. This is 
the path we follow in the next section. The intermediate 
data from Hipparcos are available in the form of one- 
dimensional 'abscissa residuals', which are the residuals 
from the st andard 5-pa rameter solutions published in the 
Catalogue (|ESA|[l997l ). By extending the 5-parameter 
model to include orbital motion, as done by the Hippar- 
cos team, these observations can be used to strengthen 
the combined solution. The nominal errors of these mea- 
surements have a median of 1.85 mas. 

4. COMBINED ORBIT AND STELLAR MASSES 

Long-baseline interferometry provides information on 
the relative orbit of the binary (with semimajor axis a) 
as well as the brightness difference, whereas the Hippar- 
cos observations refer to the motion of the center of light 
of the binary relative to the barycenter (with semimajor 
axis dphot), on an absolute frame of reference. The two 
kinds of measurements are redundant to some extent be- 
cause they both constrain the elements P, e, i, uj^, 
and T. Their simultaneous use in a global fit is there- 
fore expected to lead to a more robust solution. The use 
of the Hipparcos measurements introduces several other 
parameters that must be solved for at the same time, 
including corrections to the catalog values of the posi- 
tion of the barycenter (Aq*, AS) at the mean reference 
epoch of 1991.25, and corrections to the proper motion 
components (A/i* , Aps) and to the parallax (A^Rip) 1 . 
The formalism for modelin g the abscissa residuals follow s 
closely that describe d bv Ivan Le euweii & Evansl (|1998[ ), 
iPourbaix k. Jorissenl (|2000f ). and Uancart etldl (|2005| i. 
including the correlations between measurements from 
the two independent data reduction consor tia that pro - 
cessed the original Hipparcos observations (|ESAl[j"997f ). 
The details are reviewed in the Appendix. With the 
inclusion of the semimajor axes a and a p i 10 t there are 
13 unknowns altogether. We solve for them simultane- 
ously using stand ard non-linear least-squares techniques 
(|Press et al.lll992l p. 650). In the course of the iterations 
we noticed that one of the NPOI observations (taken on 
2004 July 21) gave an unusually large residual in both 
angular separation (+0.86 mas, or 3.2cr) and position an- 
gle ( — 1? 16, or 4. la - ). The same obse rvation stands out 
in the solution by IZavala et al. (2006) with similar resid- 

Following the practice in the Hipparcos Catalogue we define 
Aa* = Ac* cos 5 and A/x* = A/x a cos<5. 
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uals, and is one of several measurements taken at nearly 
the same orbital phase, so is not particularly critical. We 
have elected to exclude it from the solution. The relative 
weights between the NPOI and Hipparcos observations 
were assigned according to their individual errors. Since 
internal errors are not always realistic, we adjusted them 
by applying a scale factor in such a way as to achieve 
a reduced % 2 value near unity separately for each type 
of observation. This was do ne by iterations. A simi- 
lar procedure was followed bv lZavala et~aT1 (|2006l ). This 
resulted in scale factors of 0.34 and 0.24 for the NPOI 
angular separations and position angles, and 0.84 for the 
Hipparcos measurements. 

The results of this fit are given in Table Q] under the 
heading NP 01+ Hipparcos. For referenc e we include als o 
the spect roscopic solution re ported by lAikmanl (|1976|) . 
the fit by I Zavala et alj (|2006l ). and the constrained solu- 
tion by the Hi pparcos team . The position angle of the 
node given by IZavala et al.l ()2006l ) requires a change of 
quadrant to conform to the usual convention that the 
ascending node corresponds to the node at which 
the secondary is receding from the observer. The orbital 
elements from our fit are generally se en to be consis- 
tent with those of Za vala et al.l (|2006f ) (whose solution 
combines interferometry with radial velocities), but with 
considerably smaller uncertainties. 

The projection of the photocentric orbit on the plane 
of the sky along with a schematic representation the Hip- 
parcos measurements is seen in Figure HJ where the axes 
are parallel to the right ascension and declination direc- 
tions. Because these measurements are one-dimensional 
in nature, their exact location on the plane of the sky 
cannot be shown graphically. The filled circles represent 
the predicted location on the computed orbit. The dot- 
ted lines connected to each filled circle indicate the scan- 
ning direction of the Hipparcos satellite for each mea- 
surement, and show which side of the orbit the residual 
is on. The length of each dotted line represents the mag- 
nitude of the O — C residual. 2 The short line segments 
at the end of and perpendicular to the dotted lines in- 
dicate the direction along which the actual observation 
lies, although the precise location is undetermined. Oc- 
casionally more than one measurement was taken along 
the same scanning direction, in which case two or more 
short line segments appear on the same dotted lines. The 
orbit is counterclockwise (direct). The path of <f> Her on 
the plane of the sky as seen by Hipparcos is shown in 
Figure [2] The contorted pattern results from the combi- 
nation of annual proper motion (indicated by the arrow) , 
parallactic motion, and orbital motion. 

In addition to providing the parallax, the key contribu- 
tion of the Hipparcos observations regarding the masses 
of 4> Her is that they allow the measurement of the semi- 
major axis of the photocenter. This parameter can be 
expressed in terms of the relative semimajor axis (e.g., 
Ivan de Kamp| [l967') as a p hot = a(B — 13), where B is the 

2 The "O—C residuals" are not to be confused with the "abscissa 
residuals", which we refer to loosely here as Hipparcos "observa- 
tions" or "measurements" . The abscissa residuals are in fact resid- 
uals from the standard 5-parameter fit reported in the Hipparcos 
Catalogue, as stated earlier, whereas the O—C residuals (or simply 
"residuals") are the difference between the abscissa residuals and 
the computed position of the star from a model that incorporates 
orbital elements. 



mass fraction Mb /(Ma + Mb) and j3 is the light fraction 
^b/(^a + ^b)- This can be put in a more convenient form 
as 
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where q = Mb /Ma and Am is the magnitude differ- 
ence in the Hipparcos passband (H p ). Therefore, if the 
magnitude difference is known, the combination of both 
kinds of astrometric measurements allows one to solve 
for the mass ratio q. Since the total mass can also be 
determined from Kepler's Third Law given the parallax 
(TTHip), the individual masses in units of the solar mass 
follow immediately: 
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In the above expressions all angular quantities are ex- 
pressed in the same units, and the period is in units 
of side real years. The int erferometric measurement of 
Am bv lZavala et alJ (|2006l ) is therefore of crucial impor- 
tance here, and we have made use of it as an external 
constraint (along with its uncertainty) in order to solve 
for the masses. Strictly speaking, however, the value en- 
tered in the equations above must correspond to the same 
passband as the Hipparcos o bservations (H v ), w hereas 
the relevant Am reported by IZavala et al . (2006) corre- 
sponds to the passband of one of the spectral channels 
of the NPOI centered at 5500 A. We have assumed that 
the latter corresponds closely to the visual band, and we 
have applied a small correctio n to transform it to the 
Hp band using the relations by Harmanec (1998). This 
transformation depends on the individual U—B and B—V 
colors of the components, which have not been measured 
directly. But as described in they can be estimated 
using stellar evolution models (see Table [5] below) . The 
result is then AH p — 2.669 ± 0.051, where the uncer- 
tainty includes the contribution from the scatter of the 
transformation. 

The individual masses we obtain for the components 
of (j) Her are reported in Table along with other quan- 
tities derived from the orbital elements, such as the ve- 
locity amplitudes (if A , Kb) predicted for the primary 
and secondary. The uncertainties in the derived quanti- 
ties account for the correlations between the various ele- 
me nts. The tota l mass is essentially the same as derived 
bv lZavala et al.l (|2006l ). but with a significantly reduced 
error due mostly to the improved parallax. However, our 
individual masses are different than the nominal values 
they inferred, and are much closer to those expected for 
the spectral types estimated by them (see below). We 
emphasize that our mass determinations in this section 
are purely astrometric, since we have not made any use 
of the spectroscopy in our fit. Several of the orbital ele- 
ments are cons iderably more precis e in our solution than 
in the work of IZavala et alj (|2006T ). including the semi- 
major axis a that in principle depends only on the NPOI 
observations, which are common to both fits. This is 
most likely due to two factors: our rejection of one of 
the NPOI observations giving large residuals, and the 
fact that we did not use the radial-velocity observations 
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of I Aikman] (|1976l ). which appear to show some inconsis- 
tency wi th the astrometry a nd may be biasing the so- 
lution of IZavala et al.1 (|2006l ) causing larger errors. We 
investigate this issue in more detail in the next section. 
Our error in the period is also smaller despite the shorter 
time span of the observations considered in our fit (15.5 
yr) compared to theirs (40.2). We attribute this to the 
same reasons. 

5. THE RADIAL VELOCITIES OF <j> Her 

The orbital solution of IZavala et al.l ((2006) is a com- 
bined fit_oJtheir_NP 01 data with the radial velocities 
from lAikmanl (|1976f ). Although this resulted in orbital 
elements with significantly smaller errors, they instead 
ele cted to use the mass function as computed originally 
by lAikmanl (|1976f ). and with their total system mass of 
4.7 M© they inferred individual masses of 3.6 M Q for 
the primary and 1.1 M Q for the secondary. They cor- 
rectly pointed out that these values seem anomalous for 
the spectral types (which they estimated as B8V and 
A8V), the primary mass being too large and the sec- 
ondary too small. Since the secondary mass is propor- 
tional to P 1/,3 M^. 3 \/I — e 2 Kxl sini, the two main fac- 
tors that contribute to th e discrepancy appea r to be the 
larger inclination angle of IZavala et al.l ((2006) compared 
to ours, and the smaller velocity semiamplitude. Given 
that the difference in the inclination angle is not suffi- 
cient to ex plain the probl em, we conclude that the Ka 
value from lAikmanl (|1976[ ) is probably too small. 

In Fig u re [3] w e compare the original radial velocities by 
lAikmanl (|1976f l with our astrometric solution from the 
previous section (indicated with the solid line), in which 
the center-of-mass velocity has been adjusted by eye to 
fit the velocities. It is seen that the agreement is very 
good with the exception of one measurement obtained 
precisely at the descending node near phase 0.0, which 
shows a residual ^5 times larger than the scatter of the 
remaining points. This measurement happens to carry 
the highest weight for establishing the semiamplitude 
Ka, and we speculate that it may be biased and may be 
pulling the amplitude toward lower values. For ref erence 
we show also the original solution by Aik man (1976) (dot- 
ted line), in which the amplitude is lower and that mea- 
surement does not particularly stand out as an outlier. 
A new spectroscopic solution without this observation 
gives a larger semiamplitude (Ka = 2.68 ± 0.20 kms -1 ), 
as expected, although it is also a bit more uncertain be- 
cause of the lack of constraint at the maximum. We note 
that several of the other elements in this trial solution are 
also closer to our results in 21 particularly the eccentric- 
ity (e = 0.516 ± 0.036) and the longitude of periastron 
(oja = 355?6 ± 3?5), which is suggestive. 

A search for other velocities of </> Her in the litera- 
ture reveale d a set of measurements of the prim ary star 
obtained bv lAdelman. Gulliver &: Ravld (|200lD between 
1989 and 1998 that appears to be of good qual i ty, an d 
that was not used in the work of IZavala et alj (20061) . 
Examination of those data 3 shows that unfortunately 
they too lack coverage near the all-important descending 
node, and in addition there are two outliers 4 . Neverthe- 

3 W e point out that Table 1 of Adelman, Gulliver & Raylc 
(2001) containing the velocities for tj> Her has the wrong star name 
in the title: v Her instead of <j> Her. 

4 One, on HJD 2,449,134.966, has a velocity of -22.0 kms -1 that 



less, we believe them to be potential ly useful so we have 
considered them along with those of lAikmanl (|1976l ) in a 
new combined orbital solution descr ibed next. The veloc - 
ities of the secondary measured by IZavala et al.l {2006) , 
on the other hand, are of insufficient quality to be of 
use and show very poor agreement with our astrometric 
solution. 

6. NEW COMBINED SOLUTION AND REVISED MASSES 

The inclusion of radial velocities in the fit along with 
the astrometry introduces a redundancy from the fact 
that the velocity semiamplitude is constrained by both 
kinds of measurements. In addition to the 13 adjustable 
quantities considered previously, we add the center-of- 
mass velocity (7) and a velocity offset (ARV) to account 
for a possible difference in the zero points of the two ra- 
dial velocity data sets. Furthermore, since the velocities 
provide an indirect constraint on Am through eq. Q] (see 
below), we now consider Am also as an adjustable quan- 
tity instead of as a fixed constraint, and the value of 
the magnitude difference from interferometry (properly 
corrected to the H p passband) is regarded as a measure- 
ment with its correspond ing uncertainty. The scale fac- 
tor for the Aikman ( 197 6]) internal errors was fou nd to be 
1.77, and the lAdelman. Gulliver fc Ravld (|2001l) veloci- 
ties were assigned uncertainties of 0.47 kms -1 so that 
their reduced x 2 is near unity. 

The results of this 16-parameter global solution are 
listed in the last column of Table [TJ The changes com- 
pared to the previous fit are generally well within the 
errors, a sign that the spectroscopic observations we 
have added are consistent with the astrometry. The 
rms residuals for an observation of unit weight (NPOI 
p and 8, radial velocities, Hipparcos) are, respectively, 
e p = 0.13 mas, eg = 0?14, ervi = 0.25 kms -1 , 6rv2 = 
0.47 kms -1 , and enip = 1-8 mas. The individual Hippar- 
cos measurements and O — C residuals are listed in Ta- 
ble [2] In Figure [4Jd we show these residuals as a function 
of orbital phase. The O — C residuals from the original 
orbital solution by the Hipparcos team are overplotted in 
the same figure (open circles), and the two distributions 
are seen to be essentially the same. For reference, Fig- 
ure [4^, shows the Hipparcos measurements on the same 
scale before accounting for orbital motion, i.e., the ab- 
scissa residuals from the standard 5-parameter solution. 
The improvement brought about by the orbital fit is ob- 
vious. Minor systematics at the level of a few milli-arc 
seconds remain in the O—C residuals of Figure [4}z», for ex- 
ample near phase 0.7, which we believe may be reflecting 
the limitations of this data set. The residuals from the 
NPOI observations are quite similar to those reported by 
IZavala et al.l (|2006f ) and are not repeated here. The ve- 
locity residuals are given in Table [3] and Table IH along 
with the original measurements; observations that were 
rejected as described above have their residuals indicated 
in parentheses. These data are compared graphically in 
Figure[5]against the computed velocity curve from our fit. 
As an interesting test on the indirect constraint provided 
by the velocities on the magnitude difference through 
other elements, we carried out another solution in which 

is much lower than any of the others and is clearly erroneous, while 
the other, on HJD 2,448,705.595, has a velocity nearly 1.5 kms -1 
(> 3(r) higher than three other measurements at a similar phase 
that show very good interagreement. 
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the magnitude difference from interferometry is ignored. 
This fit gives AH P = 3.00 ± 0.64 mag, which is not far 
from the measured value of 2.669 ± 0.051 mag. 

7. COMPARISON WITH STELLAR EVOLUTION MODELS 

As mentioned earlier the individual masses derived 
here for the components of <p Her are in good agreement 
with those expected for the spectral types of the stars. 
A more stringent comparison may be made against cur- 
rent stellar evolution models, for example in the mass- 
luminosity (or mass-absolute magnitude) plane. The ab- 
solute magnitudes of the stars in the visual band fol- 
low from the system magnitud e (V = 4.220 ± 0.025; 
iMermilliod & Mermilliodl I1994D . the magnitude differ- 
ence determined interferometrically (AV^ = 2.57 ± 0.05), 
and our revised parallax (7r^ ip = 14.34 ± 0.35 mas), 
which corresponds to a distance of 69.7 ± 1.7 pc. Ig- 
noring extinction we obtain My = 0.100 ± 0.059 and 
My = 2.670 ± 0.074. In Figured we show the measure- 
ments against several model i sochrones from the Yonsei- 
Yale series bv lYi et al.l (l200lh ( see also iDemarque et"al"1 
20041) for ages between 100 Myr and 400 Myr and solar 
composition. The agreement is excellent, and we infer an 
age of 250 ± 100 Myr from this crude comparison. The 
assumption of solar composition here is arbitrary, and in 
fact a number of detailed analyses of the chemical com- 
positon of <p Her in the literature have generally indicated 
an iron abundance above so lar. For example, the recent 
study bv lZavala et al. ( 2006) gave a value close to [Fe/H] 
= +0.2. This, however, refers strictly to the photospheric 
composition, which is known to be peculiar in <p Her and 
other HgMn stars. In these objects some of the elements 
are enhanced by very large factors (up to 10 5 ) compared 
to the Sun, presumably due to radiatively-driv en diffu- 
sion and gravitational settling (lMichaudlll970T ). There- 
fore, the interior composition (which is the relevant quan- 
tity for the comparison with models) may be inaccessi- 
ble to the observer in this class of stars, except perhaps 
through the study of stellar oscillations, should they be 
detected. 

This apparent difficulty has motivated us to turn the 
argument around and explore the possibility of inferring 
the metallicity (as well as the evolutionary age) from the 
available observational constraints via a more systematic 
comparison with the models, as described below. In ad- 
dition to the individual masses and absolute visual mag- 
nitudes, we have considered as constraints the integrated 
colors of the system, which are among the easiest quan- 
tities to measure accur ately. The observed B — V index 
of 4> Her as reported by Mcrmilli od fc Mermilliodl (|1994f ) 
is —0.070 ± 0.005 (mean of 46 individual ground-based 
measurements) . On the other hand the Hipparcos Cata- 
logue gives the considerably redder value —0.045 ±0.003, 
based on the -Bt and Vt measurements from the Tycho 
experiment onboard the satellite along with a conver- 
sion to the Johnson system. These later measurements 
were superseded by the re-reduction that resulted in the 
Tycho-2 Catalogue (H og et aT]l2000f ). according to which 
B T = 4.155 ± 0.014 and V T = 4.216 ± 0.009. Based on 
these revised values and the conversion to t he Johnson 
system described in the Hipparcos Catalogue (ESA 1997, 
Vol. 1, Sect. 1.3) we obtain B-V = -0.055 ± 0.014. 
This is still redder than the ground-based average, al- 
though the uncertainty is perhaps more realistic so that 



the two determinations are now consistent within the er- 
rors. The weighted average, which we adopt in the fol- 
lowing, is (B-V) = -0.068 ±0.008. The U -B color of 
tfi Her as reported by |Mermilliod & Mcrmilliodj (|1994f ) is 
U-B = -0.250 ± 0.009, from the mean of 40 individual 
ground-based observations. 

To compare the six measured properties (A/a, Mb, 
My, My, U — B, B — V) against evolutionary models, 
we computed by interpolation a large grid of Yonsei- 
Yale isochrones spanning a wide range of ages (50 to 
500 Myr) and metallicities (-0.50 to +0.60 in [Fe/H]). 
Along each isochrone we interpolated the colors and mag- 
nitudes of stars in a fine grid of masses over intervals of 
±1(7 centered on our measured values of Ma and Mb- 
For all combinations of a primary and secondary mass 
taken from these two intervals we calculated the theo- 
retical integrated colors of the system, and recorded all 
combinations of the isochrone age and metallicity that 
yielded simultaneous agreement with the measured in- 
dividual absolute magnitudes and combined U — B and 
B — V indices, within their uncertainties. Figure [7] dis- 
plays all consistent models in the age /metallicity plane, 
where the point size is related to the goodness of fit as 
measured by the distance D 2 between the model and the 
observations in the six-dimensional parameter space, 

D 2 = [AM A ] 2 + [AM B ] 2 + [AM*] 2 + [AM*] 2 + 

+ [A{U-B)} 2 + [A(B-V)} 2 . 

Age and metallicity are seen to be highly correlated. The 
best match with the observations occurs near the mid- 
dle of the distribution 5 for a composition very near solar 
([Fe/H] = -0.03) and an age of 210 Myr. The indi- 
vidual masses preferred by this model are within 0.4c 
of the measured values, the agreement with the individ- 
ual absolute magnitudes is virtually perfect, and that 
with the integrated colors is ~0.1er. In Table [5] we list 
the predicted properties from this model along with the 
measured values. The individual U — B and B — V col- 
ors for the components are the ones used in <j4] to con- 
vert the interferometric magnitude difference AV into 
AH p . The best-fit isochrone is shown in Figure \§jp>. 
Based on this comparison one may conclude that the 
interior composition of 4> Her is not significantly differ- 
ent from solar, which is more or less as expected for a 
young early-type system such as this. We point out, 
however, that this relies on the colors of cf> Her being 
normal and also on the colors predicted by the stel- 
lar evolution models being realistic. Theoretical colors 
are unlikely to be in error by significant amounts, but 
we do note that there are some indications of anoma- 
lies in the spectrophotometry of HgMn stars (see, e.g., 
lAdelman|[l98llAdelman. Adelman fc Pintadoll2003l and 
references therein) of a nature similar to those seen in 
metallic-line A stars and other chemically peculiar ob- 
jects, which are thought to be connected with line blan- 
keting, particularly in the bluer spectral regions. 

8. CONCLUDING REMARKS 

5 Although a heavy element abundance very much higher than 
the Sun (even beyond the range we explored) would appear to be 
allowed by the models (but with lower significance, as seen in the 
figure), and would imply considerably younger ages, we believe 
those scenarios to be unlikely. 
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The orbital solutions in $4] and Sj6] provide an interest- 
ing illustration of the usefulness of the Hipparcos inter- 
mediate data as a valuable complement to other astro- 
metric or spectroscopic observations. The abscissa resid- 
uals from the satellite mission are publicly available for 
many thousands of binary stars distributed over the en- 
tire sky, and although a number of studies have appeared 
over the past few years that do take advantage of them, 
this is largely still an underutilized resource. In the case 
of the HgMn star <j) Her the Hipparcos measurements 
are key to establishing the individual masses of the com- 
ponents, the most fundamental of the stellar properties. 
We find good agreement between these masses and the 
magnitudes and colors of the system and current stel- 
lar evolution models, suggesting the bulk composition of 



the object is near solar. Although the faint secondary has 
now been detected, both interferometrically and spectro- 
scopically, an accurate measurement of its radial velocity 
over the orbital cycle should help considerably for reduc- 
ing the mass uncertainties, which are currently at the 8% 
and 4% level for the primary and secondary, respectively. 



We thank the referee for a number of helpful comments 
on the original manuscript. This work was partially sup- 
ported by NSF grant AST-0406183 and NASA's MASSIF 
SIM Key Project (BLF57-04). This research has made 
use of the SIMBAD database, operated at CDS, Stras- 
bourg, France, and of NASA's Astrophysics Data System 
Abstract Service. 



APPENDIX 

MAKING USE OF THE HIPPARCOS INTERMEDIATE OBSERVATIONS IN THE ORBITAL SOLUTION OF <j> Her 

The intermediate data provided with the Hipparcos catalog are the "abscissa residuals" , Au, which represent the 
difference between the satellite measurements (abscissae) along great circles and the predicted abscissae computed 
from the 5 standard astrometric parameters. The standard parameters are the position of the object (cxq, 5q) at 
the reference epoch to = 1991.25, the proper motion components (/i*, and the parallax (7THi p ). We follow here 
the notation in the Hipparcos catalog and define = aocos<5 and [i* a = /i a cos5, to include the projection factors. 
The goal of incorporating an orbital model into the analysis of the Hipparcos data is to reduce the original abscissa 
residuals to values below those obtained from the 5-parameter solution by taking into account t he orbital motion of the 
photq c enter. 6 The y 2 minimizatio n approach for doing this has been described previously by Ivan Leeuwen fc Evansl 
(|1998[ ). iPourbaix fe Jorissenl (|2000f ). Uancart et all (|2005| ). and others. We review the procedure here with additional 
details to facilitate its application in other cases. Following IPourbaix fc Jorissenl (|2000l ) the % 2 sum can be represented 
quite generally as \ 2 — 3* V~ 1 H, where 

M a 

S = Av-Y, ^A Pfc (Al) 

and 3* is the transpose of 3. In this expression Aw is the array of N abscissa residuals provided by Hipparcos, and 
dv/dpk is the array of partial derivatives of the abscissae with respect to the fcth fitted parameter. The number M of 
parameters fitted to the astrometry in our case is 12: the 5 standard Hipparcos parameters (pi = ckq, P2 — Sq, p% = /i* , 
Pi = [is, Pb = 7rmp) and 7 orbital elements (a p hot, P, e, i, wa, ^, T, represented as pk with k = 6, . . . , 12). Here 
V~ x is the inverse of t he covariance matrix of the observations, containing the abscissa uncertainties and correlation 
coefficients ()ESA| Il997. Vol. 3, eqs. 17-10 and 17-11). These are both provided in the Hipparcos catalog for each 
observation. Correlations arise because the same o riginal satellite data were reduced independently by two data 
reduction consortia (FAST and NDAC; see ESA 1997), and the results from both are typically included in all solutions. 
To be explicit, 

IVx ... \ 
v 2 



V = 



'■■ 

V ... v n J 



where each subarray Vj corresponds to a pair of FAST/NDAC measurements (F/N) and is given by 



= F (A2) 

in which a-p and ctn are the corresponding uncertainties and p is the correlation coefficient. 

The partial derivatives dv/dpk for fc = 1 to 5 in equation (|Al|) are given in the Hipparcos catalog along with the 
abscissa residuals. Th e remaining derivatives can be expressed in terms of the partial derivatives of v with respect to 
a* and S„. These are (|ESAlll997l Vol. 3, eq. 17-15) 

dv dv d£ dv dri , 
-p, — = t^^- + - = 6, . . . , 12, A3 

dpk da* a dp k dd dp k 

6 Regardless of the actual model used by the Hipparcos team to obtain the final published solution, the abscissa residuals provided are 
always those resulting from the five standard parameters as listed in the catalog. This allows complete generality in extending the model 
beyond those five parameters to include orbital motion or other motions of arbitrary complexity. 
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in which £ and 77 are the rectangular coordinates of the photocenter relative to the center of mass of the binary on the 
plane tangent to the sky at (ag, So)- The general expressions for these are 

£ = a* + fi* a (t - t ) + P Q 7r Hi p + AX 
V = S + fi s (t - t ) + P«57r H ip + AY 

in which P a and Pg are the parallactic factors. Only the last term in each of these equations is relevant in our case. 
They represent the right ascension and declination components of the orbital motion, and are conveniently expressed 
as AX = Bx + Gy and Ay = Ax + Fy. Here x and y are the rectangular coordinates in the unit orbit given by 
x = cosE — e and y = vl — e 2 sin £7, with E being the eccentric anomaly. This angle is related to the period and 
time of periastron passage through Kepler's equation, E — es'mE = 2n(t — T)/P. The symbols A, B, F, and G are 
the classical Thiele-Innes constants (see, e.g.. Ivan de Kamplll967f l. which depend only on the orbital elements a p hot, i, 
wa, and ft, and are given by 

B = — a p hot(cos lua sin f2 + sin a; a cos cosi) 
A = — a p hot(cos lua cos SI — sinwA sin Q cosi) 
G= — a p hot(^ sinwA sinf2 + cos uja cos Q cosz) 
F = — a p hot(— sin wa cos fl — cos uja sin fl cos i) 

The negative sign preceding a p hot in these equations reflects our use for this particular application of the longitude 
of periastron for the primary (wa) instead of that of the secondary (wb), for consistency with the elements of the 
spectroscopic orbit. The customary form of the Thiele-Innes constants in solving the relative orbit of a visual binary 
uses wb- Trivially the two angles differ by 180 °. 

As described bv lPourbaix fc Joriss eu (2000), the nature of the orbital solution is such that of the seven derivatives 
in equation (|A3|) the only one that needs to be considered explicitly is the derivative with respect to the semimajor 
axis, dv/dpe — dv / da p ii t- The expression for S in equation (jAlj) then reduces to 

Hphot • 



dv f dv 9£ dv drj 



_ \ - uu A 



, dpk \dcto i9a p hot dS 9a pho t 

The remaining six orbital elements (P, e, i, uja, fi, and T) do not appear explicitly but are hidden in <9£/<9a p hot = 
<9(AX)/<9a p hot and dri/da p ^ ot — d(AY) / '9a p hot- Thus, 12 parameters enter the evaluation of x 2 and are to be adjusted 
to seek its minimum: six are explicit (a p hot and Ap^, with k = 1, . . . ,5), and the remaining six are implicit. 

At each iteration towards the x 2 minimum the array of O—C residuals H is properly weighted by accounting for the 
error of each abscissa residual and the correlation between the FAST and NDAC measurements, which typically come 
in pairs. Representing one of such pairs by Ej = (rp, r N)j, it is easy to see using the definition of Vj above that the 

corresponding jth term H* 1 Sj in the x 2 sum will be 

r F \ 2 2pr F r N ( r N 



P 2 



If for some reason a measurement for only one consortium is available on a certain date (i.e., the observations are not 
paired for that particular orbit of the Hipparcos satellite), the correlation coefficient is zero and the x 2 term reduces 
to (r F /er F ) 2 or (r N /cr N ) 2 . 
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TABLE 1 
Orbital solutions for cf> Her. 



Parameter 


Aikman f 1976") 
(RVs) 


Zavala et al. ("20061 
(NPOI+RVs) 


Hipparcos 


This paper 
(NPOI+ Hipparcos) 


This paper 
(NPOI+ffipparcos+RVs) 


Adjusted quantities 












P (days) 


560.5 ± 1.7 


564.69 ± 0.13 


560.5 a 


564.783 ± 0.048 


564.834 ± 0.038 


7 (kiris -1 ) 


-16.79 ± 0.06 


-16.66 ± 0.05 






-16.642 ± 0.045 


ARV (kras" 1 ) 15 










-0.41 ± 0.12 


K A (kiris" 1 ) 


2.39 ± 0.12 


2.5 




3.02 ± 0.28 c 


2.772 ± 0.073 c 


e 


0.47 ± 0.03 


0.522 ± 0.004 


0.47 a 


0.5250 ± 0.0011 


0.52614 ± 0.00086 


uj a (deg) 


357 ± 5 


351.9 ± 2.7 


357 a 


355.0 ± 4.4 


350.8 ± 1.4 


T (HJD-2,400,000).... 


50053.7 ± 5.5 d 


50121.8 ± 1.0 


50114 ± 16 d 


50121.68 ± 0.25 


50121.43 ± 0.20 


i (deg) 




12.1 + 2.9 


Q« _l_ 1 A 
OU It 14 


q 8n -t- n 77 


o i n -t- n An 


Otiaaa fdep"^ 




189.1 ± 2.5 e 


188 ± 12 


186.2 ± 4.4 


190.4 ± 1.4 


a (mas) 




32.1 ± 0.2 




32.045 ± 0.035 


32.027 ± 0.028 


n i i- (rnas] 






9.09 ± 0.65 


8.67 ± 0.39 


8.57 ± 0.36 


Aq* (mas) 








+0.45 ± 0.32 


+0.44 ± 0.32 


Ai^ (mas) 








-0.58 ± 0.43 


-0.52 ± 0.42 


Aa/.* (mas yr 1 








-0.17 ± 0.32 


-0.15 ± 0.32 


Ai/.;; (mas vr } 








+0.01 ± 0.35 


+0.03 ± 0.34 


A7Tjji p (mas) 








+0.02 ± 0.36 


+0.07 ± 0.35 


Am (mag) 




2.57 ± 0.05 f 




2.669 ± 0.051 g 


2.672 ± 0.052 h 


Derived quantities 












K B (kins" 1 ) 




o. 1 




^ -t- o AQ 

O.DZ It U.4y 


^ 9*i 4- 9Q 


M A (M ) 




O.O 




o.U 1 It U.ZO 


Q (\K 4- n 9/1 
o.UO It U.Z4 


M B (M ) 




1.1 




1.647 ± 0.075 


1.614 ± 0.066 


q = Mb /Ma 




0.31 




0.537 ± 0.031 


0.530 ± 0.027 


M A + M B (M ) 




4.7 ± 0.6 




4.71 ± 0.37 


4.66 ± 0.34 


Ma ( mas y r_1 ) 






-25.98 ± 0.45 


-26.15 ± 0.32 


-26.13 ± 0.32 


(mas yr -1 ) 






+35.86 ± 0.48 


+35.87 ± 0.35 


+35.89 ± 0.34 


TTHip (mas) 






14.27 ± 0.52 


14.29 ± 0.36 


14.34 ± 0.35 


Other quantities pertaining to the fit 










AW 


37 


37 






36 + 18 


AWoi (p, 6) 




25 + 25 




24 + 24 


24 + 24 


AW 






76 


76 


76 


Total time span (yr) . . . 


8.4 


40.2 


3.3 


15.5 


40.2 



a Value adopted from the solution by^ikman (197^), and held fixed.* 5 Systematic radial velocity offset in the sense (Adclman, Gulliver Sz Ravlc (2001) 
annus Parameter derived from other elements in this solution, as opposed to being adjusted. 1 ^ Shifted forward by an integer number of 

cycles from the published epoch in order to match other solutions. 6 Quadrant reversed from published value/ Corresponds to a wavelength of 5500 A, and 
was derived from the intcrfcromctric visibilities separately from the orbital solution by holding the orbital elements fixed. g Converted to the Hipparcos 
passband (H p ) and held fixed in the solution, as an external constraint.* 1 Corresponds to the Hipparcos passband (H p ). 
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TABLE 2 

Hipparcos MEASUREMENTS OF <j> HER AND CORRESPONDING O — C RESIDUALS. 



HJD 






v & 





b 

V 


O- 


-C 




(2,400,000+) 


Julian Year 


(mas) 


(mas) 


(mas) 


Phase 


A 70 C A A C 1 n 

478o4.4ol0 .... 


1989 


9232 


+1.91 


2 


08 


-1 


30 


n nn a 

U.UU42 


A TQd A Q A A (\ 

4 /8o4.o44U .... 


1989 


9229 


+1.86 


1 


50 


-1 


32 


n nn a n 
U.UU4U 


478o4.o710 .... 


1989 


9238 


+3.44 


1 


29 


+0 


14 


0.0045 


/ITOi^/l 7(197 

4/8o4. /9o/ .... 


1989 


9241 


+6.07 


1 


23 


+2 


76 


n nn a q 
U.UU48 


A TCI OCT 700 

47925.7283 .... 


1990 


0910 


— 1.23 


1 


60 


+1 


74 


0.1126 


47925.0091 .... 


1990 


0905 


—3.51 


1 


55 


-0 


56 


0.1124 


a 7n /i o £; i aa. 
4/948. OlOO .... 


1990 


1536 


+2.36 


1 


62 


+0 


40 


0.1532 


4 /y45.05UU .... 


1990 


1538 


O A A 


1 


55 


+0 


48 


U. 15oo 


4 /y5o.5125 .... 


1990 


2500 


11 QQ 
— 11.00 


1 


87 


-3 


39 


U.2155 


4/y5o.o512 .... 


1990 


2501 


— 10.88 


2 


09 


-2 


95 


0.2155 


/i7nnn qq*7e 
4/999. 5o/5 .... 


1990 


2939 


— 8.12 


1 


35 


+0 


79 


0.2439 


/i Tnnn TA cro 
4/999. /452 .... 


1990 


2936 


—1.16 


1 


29 


+1 


18 


0.2437 


48Uo8.U952 .... 


1990 


3986 


+0.59 


1 


73 


-2 


24 


0.3116 


45Uo/.9/9U .... 


1990 


3983 


+1.31 


2 


01 


-1 


52 


0.3114 


45Uo8.oU88 .... 


1990 


3992 


+1.32 


3 


12 


-1 


39 


0.3120 


45Uo8.4ooo .... 


1990 


3995 


+5.24 


1 


66 


+2 


52 


0.3122 


48U55.0215 .... 


1990 


4466 


—4.62 


1 


88 


+2 


62 


0.3426 


48055.7545 .... 


1990 


4470 


—4.28 


2 


09 


+2 


95 


0.3428 


48091.5720 .... 


1990 


5450 


+16.73 


2 


30 


+4 


64 


0.4063 


48U9U.8o92 .... 


1990 


5430 


+14.25 


2 


25 


+2 


17 


U.4U5U 


4oii r .00 /o . . . . 


1990 


6156 


+4.41 


1 


69 


-1 


38 


0.4519 


<01 1 7 A 1 1 K 

48117.4115 .... 


1990 


6158 


+6.76 


1 


64 


+0 


97 


0.4520 


A 1 A d. OC7/I 

45140.20 (4 .... 


1990 


6948 


+5.40 


1 


56 


-1 


98 


0.5031 


48140. 28oo .... 


1990 


6948 


+7.56 


1 


79 


+0 


18 


0.5031 


/IOIOI 01 1 c 

48181.8115 .... 


1990 


7921 


+11.30 


1 


65 


-0 


27 


U.5ooU 


48181.8115 .... 


1990 


7921 


+9.34 


2 


13 


-2 


23 


U.50OU 


A oon. A 01 cc 
482U4.8150 .... 


1990 


8551 


—6.40 


1 


67 


+0 


62 


U.0U00 


482U4.8/08 .... 


1990 


8552 


— 1 . /9 


1 


86 


-0 


77 


U.bUo9 


/ioo/in nnon 
4824U.9U2U .... 


1990 


9539 


—2.13 


1 


63 


-0 


35 


U.o /Uo 


4824U.81oU .... 


1990 


9536 


—3.20 


1 


99 


-1 


44 


0.6705 


482o7.Uo24 .... 


1991 


0254 


— 10.06 


1 


59 


-0 


90 


0.7169 


482b0.8o8U .... 


1991 


0249 


1 n (17 
— 10.9/ 


1 


54 


-1 


80 


0.7166 


45z9o.4o4o .... 


1991 


0978 


— 10.81 


1 


88 


-1 


03 


0.7637 


/i ooriQ cm 
4829o.oU12 .... 


1991 


0982 


— 12.08 


2 


08 


-2 


30 


U. /0o9 


ion 1 c on 
48o20.152U .... 


1991 


1873 


+2.59 


1 


58 


-0 


63 


0.8216 


48o2o.Uoo8 .... 


1991 


1870 


+3.56 


1 


88 


+0 


34 


0.8214 


AQQ A A 10/IH 

48o44.184U .... 


1991 


2366 


—5.55 


1 


72 


-1 


48 


0.8535 


A OO A A AC 7n 

48o44.1o7U .... 


1991 


2366 


—4.04 


2 


17 


+0 


03 


0.8535 


48o82.U545 .... 


1991 


3403 


+4.77 


2 


41 


-0 


87 


0.9205 


45o5l.9791 .... 


1991 


3401 


+3.36 


3 


12 


-2 


28 


n < \ r\ a 

U.92U4 


45o9 ( .5598 .... 


1991 


3828 


+4.72 


1 


70 


-0 


71 


U.945U 


48o9/.5U52 .... 


1991 


3826 


+4.29 


1 


81 


-1 


14 


n a 7n 
0.94/ 9 


484o5.2521 .... 


1991 


4860 


—2.14 


1 


91 


+2 


33 


0.0147 


484o5.2451 .... 


1991 


4860 


—2.18 


2 


14 


+2 


29 


0.0147 


/IO y| C7 JflCQ 

4845 /.4U5o .... 


1991 


5466 


—2.00 


2 


09 


+2 


51 


0.0539 


/|0/| C7 /I Q 1 Q 

4545 / .4olo .... 


1991 


5467 


—3.33 


3 


00 


+1 


18 


U.U54U 


AQ AQCi OOQ7 

45459.255/ .... 


1991 


6339 


A (17 

+4.9 ( 


1 


84 


+0 


88 


0. 1104 


48489. o41 f . . . . 


1991 


6341 


+3.04 


1 


50 


-1 


05 


0.1105 


48521. 29U8 .... 


1991 


7215 


+0.32 


1 


60 


-1 


21 


0.1670 


48521. 0U8 ( . . . . 


1991 


7216 


+2.00 


1 


90 


+0 


47 


n 1 a 71 
U.lu (I 


loroo 7000 

4858o. /98o .... 


1991 


8927 


+6.38 


2 


10 


-3 


46 


n 0777 
0.2/ (1 


a cr < \ c\ r\ c\ 

4858o.90UU .... 


1991 


8930 


+9.65 


2 


18 


-0 


20 


n 0770 

0.2779 


450U (.0 (Uo .... 


1991 


9578 


—6.31 


2 


11 


+0 


34 


0.3198 


45oU (.Of ZD ... . 


1991 


9572 


—5.34 


2 


21 


+1 


30 


0.3195 


A OCH7 7<oo 

48607.7433 .... 


1991 


9582 


—7.45 


2 


17 


-0 


61 


0.3201 


Aorr\'7 707^ 

48607.7876 .... 


1991 


9584 


— 10.48 


2 


22 


-3 


64 


0.3202 


48638.0588 .... 


1992 


0412 


—0.21 


1 


78 


-1 


53 


O 070 O 

0.3738 


450O5.UUU4 .... 


1992 


0411 


—0.82 


2 


02 


-2 


15 


0.3737 


48008.1421 .... 


1992 


1236 


— 15.15 


1 


68 


-3 


91 


O. A 07fl 

U.42 /U 


48008.1951 .... 


1992 


1237 


— 13.20 


1 


75 


-1 


96 


0.4271 


A 070 C /ifiot; 

48/25. 4oo5 .... 


1992 


2805 


—2.79 


1 


84 


-2 


39 


0.5285 


A 070 c Q c 1 a 
48/25. o51o .... 


1992 


2802 


—2.23 


2 


39 


-1 


84 


0.5283 


AQT AC\ O A1 A 

45 / 4U.84 /4 .... 


1992 


3226 


— 10.82 


1 


60 


-0 


68 


0.5558 


'iOi^U.OOO'i .... 


1992 


3227 


— 1U.OU 


1 


55 


-0 


21 


u.oooo 


48798.0636 .... 


1992 


4793 


+ 1.41 


1 


80 


-1 


43 


0.6571 


48798.1008 .... 


1992 


4794 


+1.61 


1 


87 


-1 


23 


0.6571 


48832.6704 .... 


1992 


5740 


+3.87 


3 


15 


-3 


40 


0.7183 


48832.7720 .... 


1992 


5743 


+2.58 


2 


84 


-4 


68 


0.7185 


48947.7079 .... 


1992 


8890 


-5.01 


1 


95 


-0 


73 


0.9220 


48947.6419 .... 


1992 


8888 


-3.46 


3 


21 


+0 


83 


0.9219 


48947.9211 .... 


1992 


8896 


-3.18 


2 


38 


+0 


97 


0.9224 


48947.9711 .... 


1992 


8897 


+1.57 


2 


77 


+5 


72 


0.9225 


49023.4580 .... 


1993 


0964 


+9.02 


2 


36 


+3 


96 


0.0561 


49023.4832 .... 


1993 


0965 


+5.99 


4 


55 


+0 


93 


0.0561 
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TABLE 2 — Continued 



HJD 






<Jv h 


O-C 




(2,400,000+) 


Julian Year 


(mas) 


(mas) 


(mas) 


Phase 


49056.5632 .... 


1993.1870 


-4.23 


1.34 


+1.55 


0.1147 


49056.5602 .... 


1993.1870 


-4.68 


1.11 


+1.10 


0.1147 



Abscissa residuals as provided in the original Hipparcos 5-parameter solution (see text). Original uncertainties have been scaled by the factor 
0.84 (sec text). 



TABLE 3 

Radial velocity measurements of <j> Her from Airman (1976). 



HJD 

(2,400,000+) 


Julian Year 


RV 
(kms" 1 ) 


°-RV a 
(kms" 1 ) 


O-C 
(kms ) 


Phase 


38862.9364 


1965 


2784 


-13 


66 


1 


06 


+0 


39 





0676 


38910.7884 


1965 


4094 


-16 


10 





60 


-0 


04 





1523 


39328.7297 


1966 


5537 


-15 


94 





21 


+0 


14 





8922 


39637.8759 


1967 


4001 


-17 


96 





37 


-0 


11 





4396 


39657.7890 


1967 


4546 


-17 


64 





34 


+0 


27 





4748 


39658.8948 


1967 


4576 


-17 


63 





32 


+0 


28 





4768 


39897.0577 . . . . 


1968 


1097 


-16 


03 





14 


-0 


10 





8984 


39962.9980 b . . 


1968 


2902 


-13 


85 





21 


(-1 


37) 





0152 


39986.7324 


1968 


3552 


-13 


59 





27 


+0 


10 





0572 


39987.7414 


1968 


3580 


-13 


87 





11 


-0 


11 





0590 


39988.7903 


1968 


3608 


-14 


05 





18 


-0 


23 





0608 


39989.7716 


1968 


3635 


-13 


74 





19 


+0 


14 





0626 


40006.8033 


1968 


4101 


-14 


77 





23 


+0 


05 





0927 


40041.9075 


1968 


5062 


-16 


23 





21 


-0 


13 





1549 


40066.7294 


1968 


5742 


-16 


27 





30 


+0 


38 





1988 


40288.0550 


1969 


1802 


-18 


04 





67 


-0 


07 





5907 


40356.8535 


1969 


3685 


-17 


70 





11 


+0 


08 





7125 


40357.9518 


1969 


3715 


-17 


50 





21 


+0 


28 





7144 


40403.8423 


1969 


4972 


-17 


43 





16 


-0 


06 





7956 


40407.8790 


1969 


5082 


-17 


18 





53 


+0 


13 





8028 


40745.8653 


1970 


4336 


-18 


29 





51 


-0 


53 





4012 


40763.8147 


1970 


4827 


-18 


59 





50 


-0 


75 





4330 


41026.9650 


1971 


2032 


-16 


03 





51 


-0 


11 





8988 


41134.7499 .... 


1971 


4983 


-15 


36 





39 


-0 


63 





0897 


41149.7702 .... 


1971 


5394 


-14 


93 





18 


+0 


47 





1163 


41214.7197 


1971 


7172 


-17 


01 





18 


-0 


06 





2313 


41522.8695 


1972 


5609 


-17 


15 





25 


+0 


04 





7768 


41527.7210 


1972 


5742 


-17 


58 





19 


-0 


14 





7854 


41554.6734 


1972 


6480 


-16 


93 





19 


+0 


10 





8331 


41632.5999 


1972 


8613 


-13 


11 





25 


+0 


19 





9711 


41735.0464 


1973 


1418 


-16 


96 





34 


-0 


90 





1525 


41750.0084 


1973 


1828 


-16 


88 





42 


-0 


45 





1789 


41775.9408 


1973 


2538 


-17 


82 





62 


-0 


92 





2249 


41793.8875 


1973 


3029 


-17 


23 





65 


-0 


09 





2566 


41883.7190 


1973 


5489 


-18 


00 





58 


-0 


20 





4157 


41898.7083 


1973 


5899 


-17 


78 





21 


+0 


08 





4422 


41929.6650 


1973 


6746 


-18 


00 





25 


-0 


06 





4970 



a Original uncertainties have been scaled by the factor 1.77 to yield a reduced x 2 near unity in the combined fit. Observation excluded from the 
fit (see text) . 
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TABLE 4 

Radi al velocity measurements of <j> Her from 
IAdelman. Gulliver fc Rayle ( 200T|) . 



HJD 

(2,400,000+) 


Julian Year 


RV 

(kms" 1 ) 


°p,v a 
(kms -1 ) 


O-C 
(kms^ 1 ) 


Phase 


47751.7740 


1989 


6147 


-17 


00 





47 


-0 


12 





8047 


48141.7160 . . . . 


1990 


6823 


-17 


10 





47 


+0 


42 





4950 


48705.5950 b . . . 


1992 


2261 


-16 


00 





47 


(+1 


52) 





4933 


48706.9190 . . . . 


1992 


2298 


-17 


20 





47 


+0 


32 





4957 


49134.9660 b . . . 


1993 


4017 


-22 


00 





47 


(-5 


30) 





2535 


49394.9930 


1994 


1136 


-17 


80 





47 


-0 


11 





7139 


50166.9250 


1996 


2270 


-14 


20 





47 


-0 


15 





0805 


50168.0620 


1996 


2301 


-13 


30 





17 


+0 


81 





0825 


50169.0250 


1996 


2328 


-13 


50 





47 


+0 


66 





0843 


50591.7310 


1997 


3901 


-16 


80 





47 


-0 


18 





8326 


50592.9870 


1997 


3935 


-16 


50 





47 


+0 


09 





8348 


50593.9840 


1997 


3963 


-16 


80 





47 


-0 


23 





8366 


50595.9090 


1997 


4015 


-16 


70 





47 


-0 


17 





8400 


50653.9980 


1997 


5606 


-14 


80 





47 


-0 


71 





9429 


50654.9910 


1997 


5633 


-14 


70 





47 


-0 


68 





9446 


50655.8160 


1997 


5655 


-14 


00 





47 


-0 


04 





9461 


50657.0010 


1997 


5688 


-13 


90 





47 


-0 


03 





9482 


50658.9960 


1997 


5743 


-12 


90 





47 


+0 


82 





9517 


50943.8490 


1998 


3541 


-17 


30 





47 


+0 


17 





4560 


51028.8540 


1998 


5869 


-18 


10 





47 


-0 


55 





6065 



Uncertainties were not reported in the original publication, so the values listed here were derived from the requirement that the reduced % f° r 
the velocities be near unity in the combined fit. Observation excluded from the fit (see text). 



TABLE 5 

Observed and predicted properties of 6 Her. 



Object 


Mass 
(M Q ) 


My U-B 
(mag) (mag) 


B-V 

(mag) 


Observed 


Primary 

Secondary 

Combined 


3.05 ± 0.24 
1.614 ± 0.066 


0.100 ± 0.059 
2.670 ± 0.074 

-0.250 ± 0.009 


-0.068 ± 0.008 


Predicted for [Fe/H] = -0.03 and Age = 210 Myr 



Primary . . . 
Secondary . 
Combined. 



3.09 
1.590 



0.098 
2.669 



-0.267 
h0.053 
-0.249 



-0.094 
^0.271 
-0.067 



Note. — Predictions are based on Yonsei-Yalc models by Yi ct al. (2001) (sec also Dcmarquc ct al. 2004). 
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i i i i i i i i i i i i i i i i i i i i i i i i i 

10 5 0-5 -10 

Aa (mas) 

Fig. 1. — Motion of the photocenter of <j> Her relative to the center of mass of the binary (indicated by the plus sign) as seen by Hipparcos. 
The solid curve is the computed orbit from our solution combining Hipparcos and NPOI observations (Table [T}. The one-dimensional 
abscissa residuals are shown schematically with a filled circle at the predicted location, dotted lines representing the scanning direction of 
the satellite, and short perpendicular line segments indicating the undetermined location of the measurement on that line (see text). The 
length of the dotted lines represents the magnitude of the O — C residual from the computed location. One measurement with a larger 
residual was omitted for clarity. Also indicated on the plot are the location of periastron (open circle near the top) and the line of nodes 
(dashed). Motion is direct (counterclockwise). 
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100 50 -50 -100 

Aa (mas) 

Fig. 2. — Path of the center of light of <j> Her on the plane of the sky, along with the Hipparcos observations (abscissa residuals) represented 
as in Figure [T] The irregular motion is the result of the combined effects of parallax, proper motion, and orbital motion according to the 
global solution described in the text. The arrow indicates the direction and magnitude of the annual proper motion. 
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Orbital phase 



Fig. 3. — Radial velocity measurements by [Aikman (1976) compared with the predicted velocity curve from our astrometric orbital 
solution in [|4] (NPOI+ Hipparcos, solid line), in which the velo cities were not used. Phase 0.0 corresponds to periastron. The dotted line 
represents the original spectroscopic solution by Aikman ( 1 9761 V 
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Fig. 4. — (a) Abscissa residuals of <j> Her from the Hipparcos mission as a function of orbital phase. These are the residuals from the 
standard 5-parameter solution as published in the Catalogue flESAII1997D . The large scatter in excess of the internal errors is an indication 
that there is unmodeled motion in the system, (b) Filled circles represent the O — C residuals of the Hipparcos measurements from our 
orbital fit that combines astrometry and radial velocities. The scatter is now much smaller because the binary motion has been accounted 
for. The open circles (without the error bars, to avoid clutter) show the O — C residuals resulting from the orbital fit performed by the 
Hipparcos team, and are seen to be very similar to ours. 
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0.4 0.6 
Orbital Phase 



Fig. 5. — Radial velocity measurements by Aikman (1976) (filled circles) and Adclman, Gulliver & Ravle (2001) (open circles) along 
with the computed curve from our combined astrometric-spectroscopic solution. The data by |Ad'elTiam**G*uUiverfc*^l avlc ( 200j]) have been 
adjusted in this figure for the offset ARV found in our fit (see Table [T}. The center-of-mass velocity is indicated with the dotted line. 
Residuals are shown at the bottom. 
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Fig . 6. — Mass-Mv diagram for <f> Her. (a) Model isochrones from the Yonsei-Yale series bv IYi et al.l f200lTl ( see also Dcmarquc et al. 
2001) for a range of ages, as labeled. Solar composition has been assumed, (b) Best-fitting model that reproduces all six measured properties 
within their uncertainties (see text). The individual masses and absolute visual magnitudes as well as the combined U—B and B—V colors 
are simultaneously matched to better than OAcr. 
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Fig. 7. — Age and metallicity for all model isochrones by Yi ct al. (2001) (see also Dcmarque ct al. 2004) that yield a match to the 
measured properties of <f> Her (individual masses and absolute visual magnitudes, and combined U — B and B — V colors) within the 
observational errors. The size of the points is an indication of the quality of the agreement, with larger symbols representing a better fit 
to the data (see text). Solar composition is represented with a dotted line, for reference. The best fit corresponds to [Fe/H] = —0.03 and 
an age of 210 Myr. 



